Dynamics of classical wave scattering by small obstacles 
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A causality problem in the time-dependent scattering of classical waves from point scatterers 
is pointed out and analyzed. Based on an alternative model, the leading pole approximation of 
the exact scattering matrix of the square well potential, transparent expressions for the time and 
04 . position dependent Green function in a disordered medium are derived. 

(— i Wave propagation in complex media is a large and interdisciplinary field of research with many unsolved problems 

that are scientifically challenging and technologically important Jj| . Electromagnetic and acoustic waves are massless 
and described by a "classical" wave equation which is second order in the time derivative, whereas "quantum" 
particle waves are governed by the Schrodinger equation which is first order in d/dt. Scalar classical waves at a 
fixed frequency uj are thus equivalent to a particle wave at energy E ~ uj 2 . The analogies between the classical and 
quantum problems indeed lead to many cross fertilizations since solutions obtained in one field can be carried over to 
the other. For example, the discovery of localization of light in random media M was stimulated by earlier work on 
electron localization pf. Disordered systems have to be treated by statistical methods, for example by an ensemble 
average over the configurations of randomly distributed model scatterers. The (diagrammatic) perturbation theory for 
potential disorder in metals H has been helpful in understanding classical wave propagation, including (weak) wave 
localization B . An essential ingredient in these calculations is a mathematically simple yet physically meaningful and 
well-behaved model for the disorder. The point (<5-function) scatterer is often the basic building block for classical 
and quantum problems, because of its simple scattering amplitude. The associated ultraviolet divergence must be 
regularized, however Q6|5|. 

Most previous studies of classical waves in complex systems have been limited to monochromatic, steady-state wave 
£SJ ' fields, but the propagation of short pulses containing a broad band of frequencies is of considerable interest as well. 
Concrete challenges are provided by seismic waves excited by earthquakes or artificial explosions in terms of first 
arrival times and the waves which trail behind (the so called "coda" ) [|| . The power spectrum of the coda has 
received theoretical attention |9|, also in the framework of random matrix theory |f0|. Random matrix theory of 
transport, up to now mainly applied to quantum problems fllf| , identifies the "geometrical" or "universal" features of 
an observable that do not depend on the microscopic details of the system. Alternatively, however, one may pursue 
a microscopic approach trying to answer, for instance, what sort of information about the geophysical complexity 
of earth subsurface can be distilled from the total seismic trace. To this end models and approximations must be 
found which are well-behaved and sufficiently accurate over the whole complex frequency plane. In this Letter we 
point out that point defects interacting with classical waves are responsible for wrong analytical properties of the 
scattering matrix in the frequency domain equivalent to non-causality in the time domain, which appears to have 
escaped attention in the literature. We propose a different model that has all practical advantages of point scatterers 
but gives a causal response with the additional bonus that an ultraviolet cutoff is not needed. For simplicity, we limit 
attention in the following to scalar fields referring to H for a discussion of point scatterers for vector fields. After 
treating the simple one-dimensional (ID) problem, we turn to the practically more interesting three-dimensional (3D) 
space. The usefulness of the new model is illustrated by a calculation of the amplitude coda in a random medium 
with a low density of weak scatterers. 

Let us consider a classical scalar field in ID with spatially varying wave velocity c (x) . The ID problem corresponds 
to a layered system with planar sources and without lateral disorder, which is relevant as an approximation for e.g. 
the seismics of rock sediments |l2| . The wave potential 4>i(x;u)) at a given frequency w is then governed by the 
eigenvalue equation: 
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introducing the frequency dependent "scattering potential" 

V(x,u>) = ^(l-^) (2) 

which is "attractive" (V < 0) when the local velocity c (x) is smaller than the reference velocity Co (as for air bubbles 
in water) and "repulsive" (V > 0) when c (x) > Cq (as for metallic mercury in water). The physical quantity of interest 
is the retarded Green function or point source propagator (rj is a positive infinitesimal): 

G(x,x» = V a &(*;<*)&(*»• , (3) 

which is an observable for classical fields. A useful concept is the scattering matrix which connects the amplitude 
of incoming and outgoing amplitudes of a given scattering region J13J . In ID, the scattering matrix has rank 2 and 
eigenvalues Se, for the even channel £ — s and the odd channel £ = a. 

Let us consider a point model scattering potential Vg (x, w) ~ 7 (uj/cq) S (x) with a scattering strength parametrized 
by 7. A straightforward calculation gives: 

5f(c) = -^±^;5f = l, (4) 

where r,$ = 2co/7- The reflected amplitude from a pulse at time t' — and x' = L observed at time t and the same 
position x ~ L can be obtained by contour integration, collecting the pole at uj — iT$: 

G(L,L;t)= f^LJ- l (S s -S a )e^ L t c °e-^ (5) 

J ZTT 2lLU 2 

ffW>>e(t-%) 1 r r, < 

_ca e r s (t-2Z,/c o )0 ^21, _ A J 01 \ T.5 > ' U 

where is the step function. Clearly, the result is well behaved when the scattering potential is attractive (Tg < 0), 
but violates causality for a repulsive scatterer, since the reaction appears before the action. 

In order to shed light on this artifact, let us consider a finite square well insertion of thickness d and celerity (wave 
velocity) c centered at the origin and embedded into the infinite medium with celerity cq : 

VM = _ ^ _ _J _ (7) 

In the limit of vanishing d the last factor approaches the (5-function. For a finite scattering amplitude in this limit 
|l — Cq/c 2 \ must scale like d _1 , which is possible when c -C Cq. When c > co, the scattering strength necessarily 
vanishes when d — > 0. A repulsive point scatterer can only be realized by an imaginary celerity in the insertion, which 
causes the non-causality. 

Obviously, a better behaved model which retains the attractive features of the point scatterer is necessary. We 
show now that a single pole approximation of the scattering matrix [ fl3[ satisfies these requirements. The eigenvalues 
of the S'-matrix for the potential (0) reads: 

S s/a = ±e-^/- r ° (1 ~ e V 2 { ~ 0) (8) 

where Xq = (c — c ) / (c + c ) is the reflection coefficient of an isolated celerity discontinuity p4[ . The denominators 
of both S s i a vanish at 
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w„ = n— + i— log |r | = nil + iT. (9) 

All resonances (or quasi-normal modes fig]) are equally broadened by V and define the product representation of the 
scattering matrix (with the proper global phase factor) jlqj: 



6«/«(w) = -e ° j|- —2 — a (10) 
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For the attractive (repulsive) scatterer the residues vanish when n is odd (even) for t = s and n is even (odd) for 
I = a. The poles are always in the lower half of the complex frequency plane as required by causality. 

We can now introduce an approximation in which only the purely imaginary pole is taken into account, which is 
formally justified for long time scales At = t — 2L/cq ^> d/c and high reflectivities 1 — r§ <C 1. For an attractive 
scatterer 
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whereas for the repulsive scatterer the roles of the odd and even scattering channels are reversed. The time-dependent 
reflection amplitude is now well behaved: 
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Let us now turn to the 3D problem. The s-wave scattering matrix of a point scatterer at the origin V (r, to) 



7 (w/cq) S (r) reads 
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where (k*) = \J— ^k c and k c is a high momentum cut-off which is necessary to regularize the point scatterer model. 
Clearly 65 is unitary for attractive scatterers only (7 < 0) and does not have the correct analytical properties for all 
7. This do not disqualify point scatterers for studies of low velocity insertions and monochromatic illumination, but 
they are clearly unsuitable for pulsed, broad-band sources. We can find a remedy along the lines sketched above for 
ID using a spherical square-well scatterer with diameter d and velocity c. The higher-angular momentum channels 
are of higher order in Lud/co and can therefore be disregarded when At 3> d/co. The s-channel eigenvalue of the 
S- matrix read |13|| : 



S (u) = -e-' lud l CQ 1 + 
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For repulsive scatterers its poles are at u>q n = 2rnrc/d + iT with r = % log |to|. Taking into account only the lowest, 
evanescent mode 
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For attractive scatterers there is no purely imaginary pole since lo n = (2n + 1) Trc/d+iT and the lowest approximation 
involves two poles: 






The similarities and differences between ID and 3D are notable. The exact propagator in time 



AG e=0 (L,L;t) =- 

is then approximated for c > cq as: 
AG+ =0 (L,L;t)* 
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In Eqs. (O,E0) the exponential phase factors in Eqs. (ITqJTq) have been approximated by unity. In Figure 1 
we compare the exact reflection amplitudes (O) with the two approximations. The envelope of the response is 
quantitatively well represented by the lowest pole approximation. The subsequent neglect of the phase shift is a 
rather crude approximation for a single scatterer, but should be allowed in a many-scatterer configuration to be 
considered next. 

Let us now turn to the response of a disordered medium. The approximate scattering matrices derived above 
are directly related to the t-matrix eigenvalues tt — Cq (Si — 1) /2iui, which are proper partial sums in a perturbation 
theory of disordered systems with non-overlapping scatterers |)| . In the limit of a low density n of randomly distributed 
scatterers the interference between waves multiply scattered by different sites may be disregarded. When all scatterers 
are the same, the ensemble averaged Green function in wave vector <f and frequency space reads (in 3D): 
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For repulsive scatterers we employ the approximation (ttq). Transformation into time and position space gives: 
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where k = y(u>/co) — nto is the "renormalized" wave vector. In the weak scattering limit, which holds for time 
scales Ai> ITP 1 : 
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where the delayed signal or amplitude "coda" reads 



C(r;t) = - 
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This integral can be solved analytically for the first arrival 
C(c t; t) 
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The amplitude is (Lambert-Behr) exponentially damped by the imaginary part of the t-matrix. The time-delayed 
signal plotted in Fig. 2 originates from multiple scattering effects on the real part of the refractive index which 
is equivalent with a reduced effective velocity. It depends characteristically on the microscopic parameters of the 
systems, viz. the density, size, and the wave velocity (contrast) of the scatterers. Note that the long-lived incoherent 
fluctuations which arrive at the detector after diffusion (the intensity coda) contribute only to the intensity (G 2 (r; t)) , 
which is not discussed here. 

In conclusion, we have revealed the inappropriateness of point scatterers for time-dependent classical wave scatter- 
ing. An alternative model in terms of a lowest-order pole expansion is simple and well behaved. The propagator of 
a homogeneously disordered medium is calculated with the alternative model. A delayed signal or "coda" is found, 
which depends on the microscopic parameters of the scatterers. These results support speculations that the coda 
contains important information about a disordered medium that might be relevant for imaging applications. 
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FIG. 1. Time-dependence of the amplitude originating from a pulsed source and reflected by a spherical square- well potential 
in the s-channel; (a) high velocity (b) low velocity scatterer of diameter d. At = corresponds to the nominal arrival time. 
The inserts show the velocity profile in units of cq and the radius of the scatterer. Full curves denote the exact results, dashed 
curves correspond to the lowest pole approximation and in the dotted curves the phase shift has been disregarded. 

FIG. 2. Time delayed signal in the coherently propagated amplitude C(r — l,t) according to Eqs. ( |23|J24| ) in a disorderd 
medium with different scattering parameters Tr=5 (full curve), 10 (dashed), 15 (dotted), where I = cqt = V /(4-7rnco) is the 
mean free path, n the density of scatterers and T the imaginary part of the single-scatterer resonance frequency. 
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Figure 1; Bauer et al. (2001) 
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